Cooperative heterogeneous facilitation: multiple glassy states and glass-glass transition 
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The formal structure of glass singularities in the mode-coupling theory (MCT) of supercooled 
liquids dynamics is closely related to that appearing in the analysis of heterogeneous bootstrap 
percolation on Bethe lattices, random graphs and complex networks. Starting from this observa- 
tion one can build up microscopic on lattice realizations of schematic MCT based on cooperative 
facilitated spin mixtures. I discuss a microscopic implementation of the F13 schematic model in- 
cluding multiple glassy states and the glass-glass transition. Results suggest that our approach is 
flexible enough to bridge alternative theoretical descriptions of glassy matter based on the notions 
of quenched disorder and dynamic facilitation. 



Mode-coupling theory (MCT) is considered by many 
the most comprehensive first-principle approach to the 
dynamics of supercooled liquids pQ. Nevertheless, its sta- 
tus is rather problematic from a fundamental point of 
view, as the physical nature of the glass state and the 
microscopic interpretation of structural arrest are not 
yet fully elucidated. This is all the more so when we 
look at the higher-order glass singularities in structured 
and complex liquids. In this Rapid Communication, I 
show that multiple glassy states and glass-glass transi- 
tion in MCT can be understood in terms of a generalisa- 
tion of the notion of dynamic facilitation 2, 3 and boot- 
strap percolation 0J [5] . The latter is known to emerge 
in a variety of contexts including jamming of granular 
materials [B], NP-hard combinatorial optimization prob- 
lems [7] , neural and immune networks [HI [H] , and evolu- 
tionary modeling |10) . 

The formal structure of glass singularities predicted by 
MCT is encoded in the self-consistent equation 

$ = (1- $) M($), (1) 

where <1> is the asymptotic value of the correlator, and 
M is the memory kernel describing the retarded friction 
effect caused by particle caging, a physical feature asso- 
ciated with the de Gennes narrowing. We shall be con- 
cerned in the following with one-component schematic 
models in which the wavevector dependence of $ is dis- 
regarded and M is a low order polynomial. Equation ([!]) 
derived by taking the long-time of the integro-differential 
equation describing the evolution of the correlator of par- 
ticle density fluctuations-generates a hierarchy of topo- 
logically stable glass singularities, which can be classified 
in terms of bifurcations exhibited by the roots of the real 
polynomial 

Q($) =$-(l-$) M($). (2) 

Following Arnol'd notation, adopted in pQ, an Ag glass 
singularity occurs when the corresponding maximum 
root of Q has a degeneracy £, I > 2, and is defined by 



while the £th derivative is nonzero. The polynomial Q 
has always the trivial root $ = 0, corresponding to a 



liquid ergodic state, whereas nonzero values of $ corre- 
spond to a system that is unable to fully relax and hence 
can be identified with a glass nonergodic state. 

For two-parameter systems there are two basic singu- 
larities, A2 and A3, also known as fold and cusp bifur- 
cations. They have been extensively studied by using 
memory kernels given by a superposition of linear and 
nonlinear terms. In the F 12 schematic model the memory 
kernel is M($) = vi& + V2Q 2 while the F13 model is de- 
fined by M(<I>) = + vs§ 3 . The competition between 
the two terms produces a variety of nonergodic behav- 
iors: the linear term gives rise to a continuous liquid-glass 
transitions at which $ ~ e, where e is the distance from 
the critical point (e.g., e = T — T c ), while the nonlin- 
ear term induces a discontinuous liquid-glass transition, 
with the well known square-root anomaly $ — $ c ~ e 1 / 2 . 
In the F12 scenario the discontinuous line joins smoothly 
the continuous one at a tricritical point. In the F13 sce- 
nario, the discontinuous transition line terminates at an 
A 3 singularity inside the glass phase generated by the 
continuous liquid-glass transition, and therefore induc- 
ing a glass-glass transition (see Fig. 1 for a representative 
phase diagram) . The scaling form of the order parameter 
near the A3 endpoint is $ — <& c ~ e 1 / 3 , and more gener- 
ally $ — <I> C ~ e 1 " for an Ag singularity, as implied by 
the Taylor expansion of Q near the critical surface and 
Eqs. (|3]|. Thus one can observe a rich variety of noner- 
godic behaviors whose complexity is comparable to that 
of multicritical points in phase equilibria . It is a non- 
trivial result that only bifurcation singularities of type A^ 
can occur in MCT pQ. 

The F12 and F 13 scenarios were first introduced with 
the mere intention of demonstrating the existence of 
higher-order singularities and glass-glass transition, and 
then were subsequently observed in a number of experi- 
ments and numerical simulations of realistic model sys- 
tems [12-20 . It is important to emphasize that the pa- 
rameters Vi entering the memory kernel are smooth func- 
tions of the thermodynamic variables, e.g. temperature 
and density, therefore the nature of nonergodic behaviors 
predicted by MCT is purely dynamic. This is rather puz- 
zling from the statistical mechanics perspective of criti- 
cal phenomena where diverging relaxation time-scales are 
closely tied to thermodynamic singularity. It has been ar- 
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gued that this unusual situation stems from uncontrolled 
approximations. For example, the intimate connection of 
some spin-glass models with MCT has brought to the fore 
the existence of a genuine thermodynamic glass phase at 
a Kauzmann temperature T K below the putative dynamic 
glass transition predicted by MCT [21, 22 . A non-trivial 
Gibbs measure, induced by a replica-symmetry breaking, 
would therefore be actually responsible for the observed 
glassy behavior [23]. For this reason, the nature of the 
MCT has been much debated since its first appearance 
and several approaches have been attempted to clarify 
its status [2"TM30| . I will show here that the idea of dy- 
namic facilitation [3] , first introduced by Fredrickson and 
Andersen [2], offers some clues in this direction for its 
relation with bootstrap percolation provides a transpar- 
ent microscopic mechanism of structural arrest jH 131] . 
In the dynamic facilitation approach the coarse-grained 
structure of a supercooled liquid is represented by an as- 
sembly of higher/lower density mesoscopic cells. In the 
simplest version a binary spin variable, s$ = ±1, is as- 
signed to every cell i depending on its solid or liquid like 
structure and no energetic interaction among cells is as- 
sumed, H — —h^ZiSi- The crucial assumption is that 
the supercooled liquid dynamics is essentially dominated 
by the cage effect: fluctuations in the cells structure oc- 
cur if and only if there is a certain number, say /, of 
nearby liquid- like cells. / is called the facilitation pa- 
rameter and can take values in the range < / < z, 
where z is the lattice coordination: cooperative facilita- 
tion imposes / > 2, while non-cooperative dynamics only 
requires / = 1. This very schematic representation of 
the cage effect gives rise to a large variety of remarkable 
glassy behaviors, and it has long been noticed that they 
are surprisingly similar to those found in the dynamic of 
mean- field disordered systems [3"2"H3"3] . It has been re- 
cently observed that in a special case, an exact mapping 
between facilitated and disordered models with T K = 0, 
exists [35| . Since such models are so utterly different in 
their premises, it is by no means obvious that such a 
correspondence is not accidental and can be extended to 
systems with higher-order glass singularities. To clarify 
this issue, I will consider a generalization of the facili- 
tation approach |36j in which every cell i is allowed to 
have its own facilitation parameter fi (or, equivalcntly, 
an inhomogeneous local lattice connectivity) . Physically, 
this situation may arise from the coexistence of different 
lengthscales in the system, e.g., mixtures of more or less 
mobile molecules or polymers with small and large size, 
(or from a geometrically disordered environment, e.g., a 
porous matrix). In such facilitated spin mixtures the fa- 
cilitation strength can be tuned smoothly and is generally 
described by the probability distribution 

z 

where the weights {w^} controlling the facilitation 



strength satisfy the conditions 

z 

^w c = l, 0<w c <l. (5) 

C=o 

By tuning the weights one can thus explore a variety of 
different situations. Generally, one observes that when 
the fraction of spins with facilitation / = z— l,z is 
larger than that with 2 < / < z — 2, the glass transi- 
tion is continuous while in the opposite case it is discon- 
tinuous. One advantage of the facilitation approach is 
that when the lattice topology has a local tree-like struc- 
ture, one can compute exactly some key quantities, such 
as the critical temperature and the arrested part of cor- 
relation and its scaling properties near criticality. This 
can be done by exploiting the analogy with bootstrap 
percolation. Let p be the density of up spins in thermal 
equilibrium, 

for a generic spin mixture on a Bethe lattice with branch- 
ing ratio k = z — 1. As usual, one arranges the lattice as 
a tree with k branches going up from each node and one 
going down, and then proceeds downwards. In analogy 
with the heterogeneous bootstrap percolation problem, 
the probability B that a cell is in, or can be brought 
into, the liquid-like state by only rearranging the state of 
k cells above it [H ETJ EH ES] , can be cast in the form 

l-B = Bpl^ J2 (*) Bfc " <_1 (l--B)^ ,(7) 

where (•••)/ represents the average over the probability 
distribution Eq. Q. The right-hand side of Eq. Q is a 
polynomial of 1 — B, and hence the formal structure of 
Eq. Q is similar to that of schematic MCT (once 1 — B is 
formally identified with $). Singularities can therefore be 
classified according to the criteria already mentioned in 
the introduction. Nevertheless, it should be noticed that 
what would be the anolog of the MCT kernel in Eq. ^ 
can also have negative coefficients (besides containing an 
extra term of the form (1 — B) k /B), while the polyno- 
mial coefficients of the MCT memory kernel are restricted 
to non-negative ones. In fact, the sets of critical states 
which specify some Af glass-transition singularity are not 
identical to those describing the full bifurcation scenario 
of real polynomials of degree i, because the coefficients 
of the admissible polynomials Q form only a subset of 
all real coefficients. This observation means that the cor- 
respondence between MCT and the heterogeneous facil- 
itation approach is not an identity, but this still leaves 
enough room for building up models with MCT features, 
although some ingenuity may be required. It has already 
been shown, for example, that the F12 scenario is faith- 
fully reproduced in this framework J36[ [37] . To substan- 
tiate the above observation, I now will focus on the next 
higher-order glass singularity, which is the F13 scenario. 
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FIG. 1. Phase diagram for a Bethe lattice with z = 5 and 
facilitation as in Eq. |8J. The liquid-glass 1 transition is con- 
tinuous while the liquid-glass 2 and the glass 1-glass 2 are 
discontinuous. The light dashed line is the unstable branch 
of the phase diagram and shows the cuspoid structure of the 
terminal endpoint. 



Let us consider, for simplicity, a binary mixture on a 
Bethe lattice with z = 5 and 



*(fi) = C 1 - Q) s fi,2 + <l S fi,4- 



(8) 



For such a mixture, denoted here as (2,4), the probability 
B obeys the fixed-point equation: 

1 - B = p [q(l - B 4 ) + (1 - ? )(1 - B) 3 (l + W)] . (9) 

This equation is always satisfied by 1 — B = 0, while an 
additional solution with 1 — B > is found by solving 



B - 5B 2 + 3B 3 + 2qB 2 {3 - B). 



(10) 



A continuous glass transition is obtained by setting B = 1 
in the previous equation, giving: p c — 1/Aq. Using the 
relation between T and p (and setting h/k B = 1), one 
gets T c (q) = — l/ln(4(? — 1), implying that the continu- 
ous transition exists in the range 1/2 > q > 1/4. The 



discontinuous transition instead occurs when Eq. ( 10 ) is 



satisfied and its first derivative with respect to B van- 
ishes. The latter condition implies 



(95-l)(l-B) 
QB{2-B) : 



(11) 



and naturally leads to the square-root scaling near the 
discontinuous transition line. Thus the discontinuous 
transition can be graphically represented by expressing 
Eqs. ( 10 1 and ( 11 ) in parametric form in terms of B. The 
phase diagram in the plane (T, q) is shown in the Fig. [T] 
It exhibits two crossing glass transition lines, with con- 
tinuous and discontinuous nature, corresponding to a de- 
generate and generic A2 singularities. The discontinuous 
branch extends into the glass region below the continu- 
ous line up to a terminal endpoint which corresponds to 



FIG. 2. Fraction of permanently frozen spins, <&, vs temper- 
ature ratio, T/T c (q), with values of q larger than the crossing 
point. 



an A3 singularity. The location of the endpoint is found 
by simultaneously solving equation 



B = 



5 — 6g 
9- 6q' 



(12) 



(which is obtained by setting the second derivative of 
Eq. (fTrl to zero), along with Eqs.((To| and (JU|). The 



discontinuous branch located between the crossing point 
and the endpoint corresponds to a transition between two 
distinct glass states, called here glass 1 and glass 2. They 
are respectively characterized by a fractal and compact 
structure of the spanning cluster of frozen particles. The 
passage from one glass to the other can take place ei- 
ther discontinuously or without meeting any singularity, 
i.e. by circling around the endpoint (in a way much sim- 
ilar to liquid-gas transformation). The existence of two 
transitions in bootstrap percolation was first discovered 
by Fontes and Schonmann [38] in homogeneous trees and 
then found in Erdos-Renyi graphs and complex networks 
in [39j|40]. However, its relation with glass-glass transi- 
tion and MCT went unnoticed. In fact, the correspon- 
dence between Eqs. |l]) and Q naturally suggests the 
existence of further singularities in bootstrap percolation 
and cooperative facilitated models. 

Fig. [2] reports the behavior of the fraction of frozen 
spins, which is the analog of the nonergodicity parame- 
ter in the facilitation approach, when the temperature 
crosses the liquid-glass continuous transition and the 
glass-glass transition. This quantity can be exactly com- 
puted from B 36, 37J, and its expression is not reported 
here-we only notice that its general features, and in par- 
ticular the scaling properties near the critical states, are 
similar to those of B. We observe that the fraction of 
frozen spins first increases smoothly at the liquid-glass 
continuous transition and then suddenly jumps at the 
glass-glass transition. The jump decreases when q ap- 
proaches the endpoint and eventually disappears. At this 
special point, the additional condition that the second- 
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order derivative of Eq. ( 10 1 with respect to B vanishes, 



implies a cube-root scaling near the endpoint. These scal- 
ing features are exactly those expected from the F 13 sce- 
nario, and we obtain similar results for the mixtures (3,5) 
on a Bethe lattice with z = 6. 

To summarise, a close relationship exists between the 
structure of glass singularities in MCT and that of het- 
erogeneous bootstrap percolation. This allows the con- 
struction of microscopic realizations of MCT scenarios 
based on the heterogeneous cooperative facilitation ap- 
proach and provides further insights into the degree of 
universality of MCT. The role of the linear and nonlin- 
ear terms in the MCT memory kernel is played in facil- 
itated spin mixtures by the fraction of spins with facil- 
itation / = fc, k + 1 and k — 1 > / > 2, respectively. 



Their competition generates continuous and discontinu- 
ous liquid-glass transitions, while the order of singularity 
is primarily controlled by the lattice connectivity. This 
leads to multiple glassy states, glass-glass transition and 
more complex glassy behaviors. In this framework, the 
mechanism of structural arrest can be geometrically in- 
terpreted in terms of the formation of a spanning cluster 
of frozen spin having fractal or compact structure de- 
pending on the continuous or discontinuous nature of the 
glass transition. Finally, from the relation between MCT 
and mean- field disordered systems [3TJ it follows that 
quenched disorder and cooperative facilitation are two 
complementary, rather than alternative, descriptions of 
glassy matter, and this contributes to the long sought 
unifying approach to glass physics. 
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